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We introduce a novel class of higher-order, three-mode states called K-dimensional trio coherent 
states. We study their mathematical properties and prove that they form a complete set in a trun- 
cated Fock space. We also study their physical content by explicitly showing that they exhibit non- 
classical features such as oscillatory number distribution, sub-poissonian statistics, Cauchy-Schwarz 
inequality violation and phase-space quantum interferences. Finally, we propose an experimental 
scheme to realize the state with K = 2 in the quantized vibronic motion of a trapped ion. 
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1. Introduction 

Higher-order effect as well as multimode nature of quantum states play a significant role in various physical applica- 
tions. For instance, while multimode entangled states serve as necessary resources in multiuser quantum communi- 
cation network (see, e.g., 0,0,13])' higher-order character gives rise to simultaneous squeezing in several directions 
[USB or to enhancement of antibunching Q, etc. Here we deal with states that are higher-order and three-mode. 

The motivation for studying such states is twofold. Firstly, it is the academic motivation since we shall go on 
a route similar to those from 1-dimensional single-mode (two-mode) coherent states to Y-dimensional single-mode 
H, U E3, 0, E-l (two-mode 0, 0, \lTl llgj]) coherent states. Our study concerns the case of three modes 

thus generalizing the previously existing works. Our three-mode states are intrinsically nonclassical (see [l9j ] for a 
comprehensive review). It is worth mentioning here that a number of nonclassical states are intimately linked with 
a symmetry group. But this does not compulsorily hold for all nonclassical states. For example, the TV-mode sum- 
squeezed states [2(], El are always (i.e., for any N > 2) connected with the su(l,l) Lie algebra. Yet, the operators 
characterizing the Y-mode difference-squeezed state |2fl |22| forms a representation of the su(2) Lie algebra only 
for Y = 2, not for Y > 3 23]. The states we are going to study in this paper have no relation to a symmetry 
group. Secondly, the practical motivation could also be foreseen. As will be shown, a Y-dimcnsional state is in 
fact a quantum superposition of K massive distinguishable substates. Such superpositions not only help to deeper 
understand the physics world both from quantum and classical point of view, but they may also serve as resources for 
quantum computation and quantum information processing. Moreover, our three- mode states are naturally entangled 
and, as such, they potentially promise wide applications. In this connection we notice that pair coherent states (PCS) 
[241, which are of two-mode entangled character, have been found useful in testing fundamental concepts of physics 
[251 12a , 1271 |28| as well as in accomplishing quantum teleportation [23, y(| . We therefore do believe that there should 
be tasks (say, with participation of three parties like Fredkin or Toffoli gates and so on) which cannot be performed 
by single- or two-mode states but can be performed by three-mode ones such as ours. 

The states of interest are defined as the eigenstates of the operator (abc) K , with K a positive integer and a (b, c) 
the boson annihilation operator of mode a (b,c), subject to the conditions that they be also the eigenstates of the 
operators P = b^b — c^c and Q = a^a — c^c corresponding to the differences of the number of quanta in the three 
modes. It is noted that the operators a, b and c commute between each other. Denoting the state by \£,,p,q) K , the 
three following equations must be simultaneously satisfied 

(Sfe) K \t,P,q) K = £ K \t,P,q) K , (1) 
P\^P,q) K ^P\^P^) K . (2) 
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Q\€,p,q)x = <i\£,p,q)K ( 3 ) 

where £ = rexp(itp) G C and p, q are referred to as "charges" [3lL l32j] which, without loss of generality, are assumed 
to be non-negative. Equations 0) to © signify that a superposition of three-mode Fock states in the form 

oo 

\€>P>Q)k =J2 C '»\ n + < l)a\ n +P)b\ n )c ( 4 ) 

is a solution provided the coefficients c„ obey the constraint 



/ {n + q)\{n + p)\n\ K 

Cn+K -y (n + q + K)Un+p + K)l(n + K)l* { ) 

This constraint forms K independent strings {so, si, sk—i} where each of Sj consists of an infinite number of 
related coefficients Sj = {cj,Cj + K, Cj+2K, ■■■} with cj being the representative of the s 3 -string. The states q) K are 
therefore highly degenerate since the representatives Co, c\, Ck—x remain arbitrary. We now impose conditions on 
the representative Cj so that a state \£,,p,q) K j associated with the coefficients taken from the Sj-string be normalized 
to unity. (Note that the subscript Kj should not be conceived as a single index; it refers to two separate integers K 
and j.) Then, each state \£,p, q) Kj with a fixed j € [0, K — 1] has the following explicit Fock expansion 

00 cKn+j 

%p, q) Kj = N Kj (r 2 , P , q) ^ - %== \Kn + j + q) a \Kn + j + p) b \Kn + j) c (6) 

where 

p pq0 (Kn + j) = (Kn + j+p)l(Kn + j + qy.(Kn + j)\ (7) 

and 

/ °° T Km+j \ ~ 1/2 

N Kj (x,p,q) = (jc Z * ( g ) 

is the normalization coefficient. Since the overlap between states \£,,p,q) K j and |£',p', q^Kj' 1S determined by 

, c , j , , t v xxx N K3(t'*,P,q)N K j(£,,P,q) , n s 
Kj' {Z,P,q \Z,P,q) Kj = Sj? S pp , 5 qq j^-^=— (9) 

it follows that the normalized states q) Kj are orthogonal with respect to j,p, q but non-orthogonal with respect 
to £. It is worth emphasizing that for a given K there are K states q) K j (j — 0' 1) ■■■■> K — I) each of which is a 
three-mode entangled state. In other words, the eigenvalue set {£ K ,p, q} is X-degenerate corresponding to K linearly 
independent eigenstates \£,p, q) K j ■ We call them K-dimensional trio coherent states (KTCS's) because the dimension 
of the space spanned by these states is K even though each of them is embedded in a vector space of infinite dimension 
characterized by the complex variable £. In the same spirit, multidimensional entangled coherent states have also been 
investigated recently in [33j where they prove to be necessary for teleportation of quantum states of a particular kind. 
When K — 1 the KTCS reduces to the trio coherent state (TCS) introduced in [34( which is a generalization of the 
PCS 24] . So, in what follows, when KTCS's are quoted it implies K > 1. 

In the next part of this paper we first expose mathematical properties of the KTCS including proof of their 
(over)completeness. Physical properties are then studied showing that KTCS's are inherently nonclassical. Finally, 
we propose an experimental scheme to produce a KTCS in the vibrational motion of the center-of-mass of an ion 
trapped in real space by a three-dimensional harmonic potential. 



2. Mathematical properties 
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States \£,,p,q) K j with different j can be transformed from one to another by a "rotation" operator i?fn m (£,p, g) 
defined by 

s«™fcp,ff) = " Ki ^^ q \ c^«m m - i]K (io) 

N Km (r z ,p, q) 

where [x]k = x if x > and [x]k = x + K if x < 0. That is, for any l,m £ [0, if — 1], 

IC,p,9>a'; = RKim(€,p,q) \€,p,q) Km ■ (n) 

Another important property is: any KTCS can be expressed as a superposition of K phase-correlated TCS's. 
Namely, 

, e > N Kj (r 2 ,p,q) ^ / 27ri .A 

lf « Aff >« = ^(r a ,p, g ) 3 g exp J l6ci ' ,p ' 9> (12) 

where |^,p, g) = 9> 10 , N(r 2 ,p,q) = 7V 10 (r 2 ,p, g) and 

fo=£exp(^). (13) 

The expansion (|12|l can be straightforwardly verified by virtue of the identities X)m=o ex P [27T* (i — l')m/K] = KSu/ 

and X)m=o^ m l TO ) = Z^o* E^Lo^n-^+Z + 1) ■ It can be noted, due to Eq. ljl3"|l . that the TCS's superposing a 
KTCS are evenly distributed on a circle of radius r in the phase space. The inverse transformation of (|12(l is 

\ti Kj ,P,q) = N{r\p,q) £ N *\$"/ q) \^)k, • (14) 

Alternatively, making a change {£,Kj = £, exjp(2irij / K) — > x, ^ — %exp(— 2irij/K)} in (|14H and taking into account 
the identity 

\xexp(-2irij/K),p,q) Kj , = exp(-2irij j' / K) \x,P,q) K j> ( 15 ) 

which can be easily checked by use of the Fock representation © in both sides of l|15f) . we can cast l|14f) into a simpler 
but more convenient form (after changing \ back to £) as 

|C,p, g )=A(r 2 , P , g )X: 1 /^' 2 9) ^ (16) 
^ Njcj (r , p, q) 

More interestingly, it turns out that KTCS's of any two different dimensions are also related to each other. A hint to 
establish such a relation is first using (fT^f) to express a KTCS with a given K in terms of TCS's and then applying 
(I16|l to get the TCS's back in terms of K'TCS's with a different K 1 (generally K' ^ K). As a result, we arrive at 

N Kj (x,p, g) %^ exp Hfj/Q 

j'=0 j"=0 '^' ,/ 

The transformation i|17[l is most general in the sense that it recovers <|12fl when K' = l,f = and (|16|l when 
K = l,j = 0. In the special case, when K' — K, the r.h.s. of (|17[) becomes nothing else but its l.h.s., as it should be. 
In terms of the usual single-mode coherent state 

|^cxp(-M 2 /2)£^=|n), (18) 

n =0 v n - 
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the following expansion holds in general 

\^,p,q) Kj = N K3 {r 2 , Pl q) J J J ^-±N^a* /3* 7 * ,P, q) 

x exp [- (H 2 + |/3| 2 + | 7 | 2 ) /2] |a) a |/3) fc | 7 ) c (19) 

because of the completeness of the coherent state. In particular, however, a formula via three phase-correlated coherent 
states whose amplitudes a, (3 and 7 satisfy the equality a(3-f = £ seems to be more useful. That looks like this 

\Z>P,q) KJ = ^ J (r 2 ,p,g)exp[(|a| 2 + | / 3| 2 + | 7 | 2 )/2] 



f^o a h'PKj' Jo 2tt7 2vr 

x \a Kj , exp(i0)) a |/?^ exp(^ / )) i , \lK3- exp[-»(0 + 0')]) c ■ (20) 

In fact, when K = 1, j = and a = /? = 7 = ^ 1 ^ 3 , formula l|20|l reduces to formula (1) in |35|. 

We now turn to the important issue of proving that the K states q) K j (j = 0, 1, K — 1) form a complete 
set. As seen from (JSJ, any state q)kj s P ans the three-mode Fock space in which the q (p) first number states of 
mode a (b) with n a — 0, 1, q— 1 (rtj, = 0, 1, ...,p— 1) are lacking. In such a truncated Fock space the unity operator 
has to be of the form (similar matter associated with single-mode photon-added coherent states can be found in 36] ) 

00 

l P,1 = J2\ n + ^a\ n +P)b\ n )c( n \b( n +P\a( n + <l\- ( 21 ) 
n=0 

The resolution of unity of the KTCS's then amounts to existence of a weight function I 7 !/^ j(|£;| 2 ,p, q) such as to fulfill 
the following completeness condition 



A'-l „ 

E d 2 ^p,q) K3 W Kj {\e^<l)Kj{^P,q\=I P , q . (22) 

3=0 J 



To solve for W/ifj(|£| 2 ,p, q) we substitute £ = rexp(iip) and JHJ into l|22[) and integrate over <p. After the (^-integration 
the function WKj(r 2 ,p, q) is looked for in the form 

ttt ( 2 \ W{r 2 ,p,q,0) 

W Kj (r ,p,q)= 2 -, (23) 

where W{r 2 ,p, q, 0) is to be determined from the equation 

V^(a;,p, g, 0)x"dx = p pq o{n); n = 0, 1, 2, 00. (24) 



1) 



This is the classical Stieltjes power-moment problem |37|. with the set of n-th moments p pq o(n) parameterized by 
{p, q, 0}. If n in Eq. Ij24(l is extended to s — 1 where s E C, then our problem can be formulated in terms of the Mellin 
and inverse Mellin transforms |38l l3f)j which have been extensively used in the context of various kinds of generalized 
coherent states H3, HI BHH . Here, p pq0 (s - 1) is the Mellin transform, M[W(x,p, q, 0); s], of W(x,p, q, 0), 

i.e. 

j o M o(s-l)=M[W(a;,p ! g,0);s] := / W(a:,p, g, Ojaj'^da: (25) 

Jo 

and, W(x,p, q,0) in turn is the inverse Mellin transform, M~ 1 [p pq o(s — l);x], of p pq o(s — 1), i.e. 

W{x,p,q,0) =M- 1 [p pq o{s-l);x] := — / p pq0 {s - l)x~ s ds. (26) 

27TI J_ j;ocl 
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We know |44j that the solution of the simpler problem 



is given by 



W(x, l)x n dx = pi(n) = {n + 1)1, 



W{x,l) =M- 1 [pi(s-l);x] = x l e~ x . 



To make use of the known result (|28|l we twice apply the Parseval relation 
purpose here in the form 



(27) 
(28) 

which wc reformulate for our 



M- 1 [f(s-l)g(s-l);x} = J jM^fis - l);^" 1 ^ - 1); f ]• 
Now, by virtue of (|26() . I|27|) and 129|l . we have immediately 



(29) 



W(x,p,q,0) = M 1 [p pq0 {s-l):x]=M 1 [p p (s - l)p q {s - l)p (s - 1); 



dt t 

-M- X [p p {s - l)p q (s - l);t]M- l [ Po (s - 1); -] 
o t i 



dt 
T 



dr t 
-M- l [ Pp {s-l)-r]M- l [p q {s-l)--] 

T T 



M-^pois-l);- 



t q-l e -x/t 



r p-q-l e -r-t/r dT 



dt 



where in the last step we made use of Eq. (|28|) . Performing the r- integration we finally arrive at 



(30) 



W(x;p,q,0) 



t- 1 +(P+v)/ 2 e - x / t K q - p (2Vi)dt 



(31) 



where K n (t) = K^ n (t) = K\ n \{t) stands for the modified Bessel function of the second kind. There are two remarks 
to be made at this point as follows. 

• The condition IL'l'l) is that of over-completeness rather than completeness. This is due to the non-orthogonality 
of KTCS's with respect to £ and is explicitly reflected by the fact that any state \£,p,q) K j can be represented 
via the others as 



\Z,P,q) Kj =N Kj (£,P,q) J dt; \£,p,q) Kj . 



(32) 



• The solution (|31|l is not unique. According to the Carleman's (sufficient) condition |46|. our solution would be 
unique (non-unique) if the sum S = J2n=i w ith S n = [{Kn + j + q)\(Kn + j + p)\(Kn + j)!] -1 ^ 2 ™ , diverges 
(converges) . We now apply the logarithmic test 01 to examine the convergence of S. The logarithmic criterion 
says that if T = lim„_,.oo ln(5 n )/ln(n) > —1 (< —1) then S diverges (converges). To calculate T we rewrite S n 
in terms of Gamma functions, S n — [T(Kn + j + q + l)T(Kn + j + p + l)T{Kn + j + l)] -1 / 2 ™ , and then use 
the Stirling's formula, T(Kn + I) ~ y/2n exp(— Kn)(Kn) Kn , which is valid for large n and finite I (I <C n). As 
a result, we obtain lim^a, Ir^S^) = — 3A"ln(n)/2 and, hence, T = —3K/2 < —1 for any positive integers K. 
This proves the non- uniqueness of our solution 1311) . 

To end this section we remind from that the KTCS's are normalized by the normalization coefficient 
NKj(r 2 ,p, q) determined by JHJ. The states are also continuous in the variable £ because \£,p, q) Kj — |£',p, q) 



1 - Re 



(kj (C,p,q \t,p,q) 



'Kj 



Kj 



tends to vanish when |£ — £'| — > 0, as is evident from (JSJ and 10. Furthermore, 

as followed from l|23|) and l|31|l . the weight function WKj(\^\ 2 ,P,q) in 11' I'D is guaranteed to be positive (thanks to 
positivity of the Bessel function K n (t)). Therefore, the minimal set of conditions |48ll49| required for an ensemble of 
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TCS 



FIG. 1: 3D bar plot showing oscillation of the number distribution, P n (£,p,q, K,j), for r — 30, p = q = and different 
combinations of Kj. For comparison, also plotted is the case of TCS which does not oscillate at all. 



states to be coherent is met for the KTCS's. That explains why "coherent" arises in naming the states \^,p,q) K j as 
"K- dimensional trio coherent states" of which "coherent" must be family name and "K- dimensional trio" first name. 

3. Physical content 

In this section we explore the physical content of the KTCS's. The probability of finding I (m,n) quanta of mode a 
(b, c) in state \£,,p, a )xj IS gi ven by 



with 



Pimn(£,,P,q,K,j) = \ Kj (£,p,q\ l,m,n) abc \ 2 = P n (^,p,q,K,j)S m<n+ pS l>n+g (33) 
P n (£, p, q, K, j) = K > K ,Pq > K K ' (34) 

PpqO[n) 

where I(x) = 1 if x is an integer and I(x) = if x is a non-integer. While the Kronecker symbols in l|33|l reveal the 
entanglement of the modes in state \£,,p,q) K j , the function in (|34|l indicates that for any K > 1 the number 

distribution of state \£,p, q) Kj sun?ers an oscillation as shown in Fig. 1 in which the TCS is also plotted that does not 
oscillate at all. 

Further information on inherent quantum statistics of state \£,p, q) K a can be obtained from the Mandel parameter 
M x for mode x 

(n x 2) ) - (n x ) 2 

M x = ± (35) 

(%) 

where with x the annihilation operator of mode x and the expectation value of the factorial moment of the 

number operator (ri2^ = (jlmLoi^x — ni)^ is derived for our KTCS's in the form 
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FIG. 2: Mandel parameters M a = M b — M c versus z — r 2 for j — p = q — and different _R" indicated near the curve. For 
comparison, also plotted is the case of TCS (dashed curve) which remains sub-poissonian in the whole range of z. 




FIG. 3: Mandel parameter M b versus z for j — q — while (a) p = 1 and K = 2, 3, 4; (b) JiT = 4 and p = 1, 2, 3, 4. The dashed 
curve in (a) is the case of TCS which exhibits no maxima. 



(«?») = ^£(4). <*» 

where A'xj = Nxj(z,p, q) and z = r 2 . An evident factor that makes difference between modes is the charges p, q. For 
p = q = all the three modes are "identical" . For p > q = (q > p = 0) mode a (6) and mode c behave identically 
which are however distinct from mode b (a). For p = q > modes a and 6 are similar while mode c is dissimilar. Only 
when p ^ q and each of them acquires a positive integer the three modes are well distinguishable one from another. 
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20 40 60 



z 

FIG. 4: Mandel parameters (a) M c and (b) Mb versus z for K = 3, j = 0, p = 2 while g = and q = 2. 




FIG. 5: Mandel parameters M a , Mb and M c versus z for K = 4, j = 0, p = 1 and g = 2. 



Use of (|3*H)) in (|3"5|) yields explicitly 



2z 2 N Kl N 



N' 2 



qN 2 



2zN Kj N' Kj - q N 2 Kj 



(39) 



2z 2 N Kn N 



M h 



pN 2 



Kj 



2zN Kj N' Kj 



(40) 
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FIG. 6: Mandel parameters (a) M b and (b) M c versus z for K = 4, j = 0, 1, 2, 3, p = 1 and q = 2. 



* NkjN 



N 



Kj) 



(41) 



where iV^ = dN Kj /dz and iVj^. = d 2 N K] /dz 2 . 

In Fig. 2 we plot M c as a function of z for j = p — q = and different if. Contrary to the TCS for which mode c 
remains sub-poissonian (i.e. M c < 0) in the whole range of z, for KTCS's the mode is super-poissonian (i.e. M c > 0) 
at small values of z but then becomes sub-poissonian when z increases. The crossover point z cross at which the 
statistics changes from super- to sub-poissonian moves to the right for greater values of K. At large values of z the 
mode gets more antibunched in a higher dimension. 

For p > q — the shape of M c = M a looks like that in Fig. 2 but for a fixed K the crossover point z cross shifts 
to the large- z side with increasing p. For instance, when K = 3 numerical calculations give z cross — 7.5628, 12.0114, 
16.3108 and 20.5606 for p = 0, 1, 2 and 3, respectively. Concerning mode b, its behavior is qualitatively different 
as displayed in Fig. 3. Contrary to the TCS for which Mb monotonically increases with z, for KTCS's it exhibits a 
maximum which may be located in the super-poissonian domain if the dimension K is high enough (e.g., if K > 4 
when p — 1 as illustrated in Fig. 3 (a)). The effect of the charge p is to pull down and right-shift the whole curve 
(see Fig. 3 (b)) so that for a given K the mode can be made entirely sub-poissonian (if it was not so) by setting p 
large enough (see, e.g., Fig. 3 (b) for K = 4 : Mb < in the whole range of z when p > 3). 

For p = q > 0, as mentioned above, M a = Mb ^ M c . Compared to the situation p > q = the following properties 
hold. For mode c we find the relation M c (K,p = q > 0) > M c (K,p > q = 0) in the whole range of z except for 
z = at which M c (K,p = q > 0) = M c (K,p > q = 0), as seen from Fig. 4 (a). For mode b (a) we find that 
m&x[Mb(K,p — q > 0)] — maxf-M^-J^p > q = 0)] but the maximum of Mb(K,p > q — 0) appears "earlier" than that 
of Mb(K,p = q > 0) when z is increasing (see Fig. 4 (b)). 

For p > 0, q > and p ^ q each mode develops its own dependence on the parameters and are well distinguishable 
as Fig. 5 shows. 

The specific feature of KTCS's is their degeneracy as explained in Section 1. A given manifold characterized by a 
fixed K > 1 consists of K eigenstates enumerated by the parameter j = 0, 1 , 2, K — 1. States with the same K 
but different j differ in their number distribution which is dictated by the Kronecker symbols in Eq. I|33|l and by the 
function J(^j^) in Eq. (|34|) . These mean that one cannot find n quanta of mode c (b, a) in state \£,p, q) Krj if n — j 
(n — j — p, n — j — q) is not a multiple of K (see Fig. 1, for verification). States with the same K but different j 
differ also in their quantum statistics. So far (in Figs. 2 to 5) we have treated only j = in which case, for arbitrary 
K, p and q, we have M c (j = 0, z ^ 0) ^ K — 1 and M a b(j — 0, z — > 0) — > —1. For j > 0, however, all the three 
modes are highly antibunched at small values of z independent of K, p and q, i.e. M a ^, c (j > 0, z — > 0) — > —1. Also, 
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at small z a state with greater j has higher degree of antibunching, i.e. its number distribution profile is narrower. 
Those observations are demonstrated in Fig. 6. 

Another intriguing figure of merit is strong correlations between modes of KTCS's. Such correlations are expected 
to be nonclassical. To verify it we examine the Cauchy-Schwarz inequality (CSI) [5l| which describes a classical 
correlation between two modes x and y 



J xy = (n (2) ) (nW) - (n x n y ) 2 > 0. 



(42) 



General expressions for expectation values of products of two factorial moments n x ' can also be analytically 
derived for our KTCS's in terms of the normalization coefficient ijSJ. As a result of derivation, we arrive at 



jrr, 

-p/v 2 



l+p-q_ 



dz l 



N 2 

, Kj , 



(43) 



''N'- 



Kj 



dz r 



dz l 



N h, 



(44) 



l(0ft(m)\ = J-V N l d_ 



Use of and - gSJ in g2J) yields explicitly 



dz< 



N h 



(45) 



Jab 



N K 3 3 {pq(l -p- q)N 3 Kj + 24z 3 A^- 

~2z 2 N Kj N' Kj [(2 + 7(p + q) - (p - qf) N' Kj + izNfa] 

+ 2zN 2 Kj [6pqN' Kj + z (p + q - (p - q) 2 ) N'^] } , 



(46) 



Ja 



2z 2 N~ 3 . {!2ziV; 3 



Kj 



9(1 



Q)N 2 Kj 



N" 
ly Kj 



NK 3 N' Kj [(2 + 7q-q 2 )N' Kj + 



* zN Kj. 



(47) 



J bc = 2z 2 N K 3 {l2zN^+p(l - P )N 2 Kj N^ - N Kj N' K] [(2 + 7p - p 2 ) N' Kj + izN'^] } 



(48) 



It is well-known that a nonclassical correlation violates the CSI |52|. i.e. it makes J xy < 0. To assess degree of CSI 
violation we scale J xy to (n x n y ) 2 , i.e. we use a quantity G xy defined by 



G xy — 



>xy 



{n x n y ) 



(49) 



as a measure of CSI violation. For states with j > (i.e. j — 1, 2, K — 1) we find that the CSI is always violated. 
However, for j = 0, though the TCS also always violates the CSI, the KTCS's do not. The violation depends on the 
type of correlation and the value of charges p, q. The simulation reveals that G a b is always negative but G ac and Gbc 
may be positive depending on the charges. More concretely, when both p and q are not greater than one we find that 
G xy < Vie, y, K, z, i.e. the CSI is fully violated. However, when p (q) = 1 and q (p) > 2 the quantity Gbc (G ac ) 
remains always negative but G ac (Gbc) becomes positive at small values of z and the higher the dimension K the 
wider the z-domain within which the CSI is not violated, i.e. there is a partial violation of the CSI. Those behaviors 
are illustrated in Fig. 7. Finally, when both p > 2 and q > 2, both the quantities Gbc and G ac violate the CSI only 
partially in the sense mentioned above. 

We next study the phase-space characteristics of the KTCS's. For that purpose we consider the three-mode Q- 
function defined as 



Qfp q {a,P,l) = — abc{a 1 P 1 "f\tp,q) Kj 



(50) 
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FIG. 7: Correlation measures (a) Gj c and (b) G ac versus z for j = 0, p = 1, q = 2 and different if as indicated near the curve. 
For comparison, also plotted is the case of TCS (dashed curve) which remains negative in the whole range of z. 
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FIG. 10: Contour plots of Q = Tv 3 QfJ q (x, y) for £ = 12, p = g = 0, if = 3 with j = ( left), 1 (middle) and 2 (right). 



where a,/?, 7 £ C and |a,/3,7) a bc = l a )a|/3)fc|7) c with \a) a , \/3)b, \l)c the usual coherent states (fTHf) of modes a, b, c. 
This function is non-negative definite, bounded and normalized to unity 



/ / J QUk*>P>i)* a *P*7 = 1 ' ( 51 ) 

Generally there are six variables associated with the real and imaginary parts of a, (3 and 7. For visualization let us 
confine ourselves to a subspace determined by a = j3 — 7. In that subspace the Q-function for KTCS's is calculated 
to be 



S^'V) = -N 2 KJ (^P^) e-^+^ + y 2 )^ N^(x-iy) 3 tP ,q) 



(52) 



with x = Re(a) and y = Im(a). We represent in Figs. 8 and 9 the function QfJ q {x,y) and its contour plot for 
K = 2, j = and 1. The figures clearly manifest signatures of Schrodinger-cat-like states: constructive (Fig. 8) and 
destructive (Fig. 9) interference fringes between bell- like peaks (there are 3K bells for a given K) . The K = 2 state 
with j = (j = 1) is called even (odd) trio coherent state We also present in Fig. 10 contour plots of the 

case of K = 3 with j = 0, 1 and 2 which shows nine interfering bells. Transparently the pronounced interference 
fringe structures are j-dependent. It is the phase-space quantum interferences between different TCS's (instead of 
being simply added) that bring about copious nonclassical signatures of the KTCS. 



4. Physical realization 

After having studied physical properties of KTCS's we proceed to find ways to realize them. In this section we are 
concerned with the context of ion trap. Since ions can be trapped very efficiently and their entanglement with the 
environment is extremely weak, tra ppe d ions have advantages for many purposes such as preparing various types 
of nonclassical states (see e. g . |5 6l l57l l58l l59l l60l l6ll 162^). simulating nonlinear interactions j6j|, demonstrating 
quantum phase transitions |64tl65|. establishing quantum search algorithms 66] and so on. The most promising merit 
of trapped ion systems is perhaps the possibility to implement scalable quantum computers [67} in which a number 
of ions are involved [6^, |63, [Til • Nevertheless, many tasks can still be done even with a single ion. For instance, a 
controlled-NOT quantum logic gate can be performed just by a single trapped ion 0, IT^, El El : the target qubit 
is stored in the ion internal electronic states while the external quantized motional states, i.e. the phonon states, of 
the same ion serves as the control qubit. Here we propose an experimental scheme to generate KTCS's with K = 2 
in the vibronic motion of an ion which is trapped in real three-dimensional (3D) space. The situations corresponding 
to ID and 2D were already considered in [73 and for the single- and two-mode case, respectively. 

We first trap a two-level ion of mass M by a 31? isotropic harmonic potential characterized by the trap (phonon) 
frequency v. Let a, b and cbe the phonon annihilation operator in the x-, y- and z-axis, respectively. The ion is then 
simultaneously driven by fourteen traveling- wave lasers (compare with the TCS case 1351 ) in the resolved sideband 
regime. The first thirteen lasers are all tuned to be resonant with the sixth red motional sideband, i.e. their frequency 
is to — A — 6v with A the energy difference between the two electronic levels of the ion. For our purpose it is 
essential that the lasers be judiciously configured. Namely, we shine the 1st (2nd, 3rd, 4th, 5th, 6th, 7th, 8th, 9th, 
10th, 11th, 12th and 13th ) laser along the direction connecting the coordinate origin {x,y,z} = {0,0,0} to a point 
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{x,y,z} = {1,1,1} ({1,-1,1}, {1,1,-1}, {1,-1,-1}, {1,1,0}, {1,-1,0}, {1,0,1}, {1,0,-1}, {0,1,1}, {0,1,-1}, 
{1,0,0}, {0,1,0} and {0,0,1}). As for the 14th laser, it must be in resonance with the electronic transition, i.e. 
its frequency is equal to A, but its propagation direction is unimportant. The Hamiltonian of the ion-phonon-laser 
system is {%= 1) 



where 



H = y<T 2 + v{a)a + b f b + c+c) + H r , 
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H m t = ^2 cxp[-i(ujit + ipi) + ik / R/]er + + H.c. 

i=i 



(53) 



(54) 



with Qi the Rabi frequencies, wi,2,...,13 = to, u)i& — A, (pi the phases, R; the position operator along the laser direction 
determined by the wave vector k;, a z — |e) (e| — \g) (g\ , er+ = al_ = \e) (g\ and \g) (|e)) the ion's electronic ground 
(excited) state. Assuming ki — k VI there is a single Lamb-Dicke parameter r\ = kj\j2Mv which measures the 
localization of the ion relative to the laser wavelength. In terms of r\ we obtain, in an interaction picture with respect 
to Hq = H — Hi n t, the interaction Hamiltonian of the form 



Hint — c 



00 I ~2\m 



(-v 2 y 



13 

2=1 



(m + 6)! 



'14 V / 



4'' 



a+ + H.c. 



In the Lamb-Dicke limit 77 <C 1 we can retain only the m = term in l|55(l and reduce it to 



(55) 



H 



int 



6 13 

'■ 2=1 



o-+ + H.c. 



Thanks to our configuration for the lasers the operators Ai,2,...,i3 are expressed through a, & and c as 

A\ o = a ± 6 + c, A3 4 = a ± & — c, 



(56) 



(57) 



A 5i6 = a±b, A 7 a = a±c, A 9iW = b±c 



(58) 



An = a, A12 = b, A i3 = c. 
If we control the laser intensities and phases in such a way that 



(59) 



^1,2,3,4 — 2^5,6,7,8,9,10 — 2^11,12,13 — ^, 



^1,2,3,4 = ^11, 12,13 = <P + 1", ^5,6,7,8,9,10 = V, <Pli = 7T, 

then use of Eqs. (O - lO in lO yields 
thanks to the identity 



cr + + H.c, 
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Y, Af - 2 £ At + 4 A n = 360(S>e) 2 



(60) 
(61) 

(62) 
(63) 
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FIG. 11: 3D bar plots of phonon number distribution Il n at different values of Ft for £ = 8, p — q — 0, £ = 0.02 and I = 3. (a) 
w = 0.5: no KTCS's are generated since in the long-time limit II n is non-zero at both even and odd n, i.e. it does not display 
required oscillation, (b) w = 0: generation of state \£,p,q) 20 is justified in the long-time limit by coincidence of II n , Eq. I72L 
with P n = P n (£, p, q, 2, 0), Eq. (1341 . (c) w = 1: generation of state \£,p,q) 21 is justified in the long-time limit by coincidence 
of Tin, Eq. with Pn — Pn 



The Hamiltonian (|62l) is the central result for the physical realization under consideration. The new parameters £, £ 
appearing in (|62|) are controllable and given simply by 

C=-£-exp(-i<p), (64) 



e 2 = -^expM. (65) 

Since the trapped ion is well isolated the dominant channel of decoherence is via the spontaneous decay. Then the 
time evolution of the system density operator p is governed by the following master equation 

d P ;ru l i r ( < 7 + <J -P P< J + (J -\ (( . R ^ 

— = -i[H int ,p\ + T [a-p<J + — J (66) 

where T accounts for the spontaneous decay rate of the ion being in its electronic excited state |e). 

It is clear that the system ceases to evolve when the ion's fluorescence stops. Hence, the "dark" stationary solution 
of Eq. I|66|) in the long-time limit has the ansatz 



p(oo) = i 5 ) m xyz m (si (e?) 

where l*)^^ solely determines the phonon state. Setting dp(oo)/dt = in Eq. (|FJFJ|l and using properties of the 
operators a± we arrive at an equation for \^) xyz , 



(^ 2 m xyz = ew xyz , (68) 

where £ 2 depends on the experiment parameters via (|65p . 

For \^} xyz to be a desired KTCS \£,p, q) 2 Q or \£,p, q) 21 we need to prepare an appropriate initial state. Furthermore, 
the generation time is influenced by all the parameters involved. To see these let us simulate Eq. (|66|) by the Monte 
Carlo Wave- Function approach |72| with an initial state of the form 



|<M0)) = \e) (V^\Vio) xyz + VT^ \*n) xyz ) (69) 

where < w < 1 and 



l*«>»v* = \ 21 + 1 + J')x I2Z + p + j) y \2l + j) z 



(70) 
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rt rt 

FIG. 12: Fidelities (a) Fo and (b) Fi as a function of Ft for different sets of parameters: Curves 1, 2, 3 and 4 correspond to 
p = Q,l} = {10,0.005,2,0}, {10,0.005,0,0}, {6,0.005,0,0} and {10,0.002,0,0}, respectively. 

with a non- negative integer I and j = or 1. Since the dynamics under consideration does not mix the parity, at time 
t > we have 



IM*)) = E ( G ™(*) Iff) + l e » (V^I^O)^ + VT^ \*ml) xyz ) (71) 

m=0 

where G m and _E m are some time-dependent coefficients to be simulated. In Fig. 11 we plot the phonon number 
distribution 



U n (t) = ($ w (t)| n + q) x \n + p)\n) i 



(72) 



at different times. Obviously, for < w < 1, n n may be nonzero at any n. In particular, for w = 0.5 there is no 
oscillation at all in II„ in the long-time limit (see Fig. 11 (a)). Thus, a desired KTCS \£,p,q) 20 or \^,p, q) 2l results 
only when either w = or w = 1 as seen from Fig. 11 (b) or Fig. 11 (c). These figures show that in the course 
of time the phonon number distribution H n (t) changes gradually and finally tends to coincide with P n (£,p, q, 2, 0) if 
w = (Fig. 11 (b)) or with P n (£,p, q, 2, 1) if w = 1 (Fig. 11 (c)), implying successful generation of the target state 

lf,p,g) 2 o or 1^9)21 • 

To assess quality of the resulting state as well as the generation time of our scheme we examine time-dependence 
of the fidelity 



^■W = I<^WU,P,?> 2J | 2 - (73) 

In Fig. 12 we represent Fj as a function of Tt for various sets of parameters. Although initially the transient behaviors 
happen differently for different parameters' sets, after some period of time a stationary regime is inevitably established 
in all cases. It also follows from the figure that the generation time, i.e. the time needed to reach the stationary 
regime, is shorter for greater p, q (compare Curves 1 and 2), f (compare Curves 2 and 3) as well as ( (compare 
Curves 2 and 4). For a wide range of parameters we find out that 1 — Fj becomes less than 10~ 3 for Ft > 200. The 
desired KTCS is thus produced with high purity within a relatively short period of time (compared, e.g., with the 
TCS generation time |35(). The KTCS obtained in our scheme is also stable since its appearance is accompanied by 
the "dark" ion which is to be found in the ground state \g) and therefore does not interact with the laser fields any 
longer. 

5. Conclusion 

We have introduced and studied KTCS's which are both higher-order and multimode. Being superposed of phase- 
correlated TCS's, a KTCS is indeed a new physical state as dictated by the quantum mechanics superposition 
principle. The novel nonclassical features of KTCS's, as compared to TCS's, were demonstrated in detail and an 
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experimental generation scheme for K = 2 was also presented in the ion trap context. Since these KTCS's are entangled 
states, they do promise potential implementations in quantum information processing and quantum computation. 
In particular, their three-mode nature would be crucial for tasks such as quantum controlled teleportation or/and 
quantum telecloning of certain types of continuous-variable states as well as inequality-free "one-shot" tests of local 
hidden- variable theories, etc. Such applications of KTCS's are being under way and will be reported elsewhere. 
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